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We investigate the role played by nonequilibrium dynamical screening in the thermal- 
\ ization of carriers in a simplified two-component two-band model of a semiconductor. The 

. main feature of our approach is the theoretically sound treatment of collisions. We abandon 

' Fermi's Golden rule in favor of the Schwinger-Bakshi-Mahantappa-Keldysh (nonequilibrium 

field theoretic) formalism as the former is applicable only in the long-time regime. We also 
» I ■ introduce the concept of nonequilibrium dynamical screening. The dephasing of excitonic 

f> f quantum beats as a result of carrier-carrier scattering is brought out. At low densities it is 

found that the dephasing times due to carrier-carrier scattering is in picoseconds and not 
femtoseconds, in agreement with experiments. The polarization dephasing rates are com- 
puted as a function of the excited carrier density and it is found that the dephasing rate 
for carrier-carrier scattering is proportional to the carrier density at ultralow densities. The 
scaling relation is sublinear at higher densities, which enables a comparison with experiment. 



< 



> 

Q\ • I. INTRODUCTION 



O 



The relaxation of hot electron distributions via emission of LO optical phonons has been studied by 
a number of authors. However, the role of carrier-carrier scattering has received less attention for two 
main reasons. The theory of screening involved in the proper description of carrier-carrier scattering is 
poorly understood. Any reasonable attempt to include carrier-carrier scattering is computationally intensive. 
Secondly, the perception that carrier-carrier scattering is significant only for large carrier densities has 
resulted in few attempts at accurately modeling the phenomenon. In particular, authors in Refs. have 
conducted theoretical studies of the relaxation of hot electrons when the emission of LO optical phonons 
is the dominant mechanism. The authors in Ref. [5]-Q have attempted to investigate the role of carrier- 
carrier scattering in momentum and energy relaxation. On the experimental side, authors in Refs. p|-|l5| 
have investigated the relaxation of hot carriers in bulk GaAs as well as quantum well structures such as 
GaAs/AlGaAs. A number of experiments have been performed in the recent years on the density dependence 
of the scattering rate of free carriers which provide information about the role of carrier-carrier scattering. 
Prominent among them are the experiments by authors in Refs. |5| and 16 2^. The use of nonequlibrium 
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' Green function techniques in the study of relaxation phenomena is not widespread. However, attempts have 

been made in the recent past. Notable among those are the study of memory effects due to relaxation by 
phonons by authors in Ref. |35| and an attempt by authors in Ref. [36] to model the nonequlibrium aspects of 
screening of the Coulomb interaction, similar to the approach outlined here. 

It has been pointed out in Ref. |^ that changes in the scattering rates with carrier density are solely due 
to coulomb scattering of carriers. Therefore an investigation of the density dependence of carrier-carrier 
scattering should provide a good testing ground for theories of coulomb scattering. The authors in Ref. ||-|7] 
have developed a theory of carrier-carrier scattering using the quasi-classical Boltzmann equation coupled 
with Fermi's Golden rule for evaluating the the collision terms. Objections to the use of Fermi's Golden rule 
may be raised, however. Firstly, Fermi's Golden rule is derived from time-dependent perturbation theory as 
an asymptotic time approximation which translates into a_strict energy conservation rule. Quantitatively, 
the time resolution of a typical photon echo measurementEEl is of the order of At ~ 10/s.. This translates 
to an energy uncertainty of an order of AE ~ SOmeV. This is substantial when compared to the plasmon 
energy at densities of the order of po ~ 2 x 10 16 /cm 3 or the LO phonon energy which is ujlo = 36meV. 
At these densities the plasmon energy is lj p ~ 6.2meV. We therefore abandon the familiar rule in favor of 
a more elaborate explicitly nonequilibrium field-theoretic formulation which we shall discuss subsequently. 
There is at least one other significant conceptual drawback in the present theories of Coulomb scattering in 
carrier thermalization. We are primarily interested in the dynamics that occur within 10 ps. of the switching 
on of the external field. Within this time scale, the system is in a highly nonequilibrium state. Therefore it 
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does not, in particular, make sense to use the naive generalisation of the equilibrium dynamically screened 
coulomb interaction which essentially involves replacing in the formula for the RPA dielectric function, the 
equilibrium carrier densities by nonequilibrium ones. This is because when time translation invariance is 
broken as is the case here, the RPA (random phase phase approximation) as it is commonly understood, is 
no longer valid as it is not possible to transform to the frequency domain. ■Qux.field- theoretic formulation 
involves the use of Schwinger-Bakshi-Mahantappa-Keldysh-Kadanoff-BaymOLZI approach in deriving the 
collision terms (henceforth the phrase 'collision terms' will refer to the right side of the Boltzmann equation. 
A good introduction to this approach may be found in Ref. |7|.) We assume that the collision terms in the 
Boltzmann equation can be expanded in a power series of the external (time-dependent) vector potential. 
This approach is simliar to linear (or quadratic) response theory. The main difference being that we solve 
the full Boltzmann equation numerically after the collision terms have been obtained so that the densities 
and currents no longer depend only on finite powers of the external fields. We also obtain explicit formulas 
for the screened coulomb interaction in terms of the external fields. These formulas clearly indicate the 
nonequilibrium nature of the dynamical screening. This screened potential also produces a finite collision 
cross-section in both two and three space dimensions. 



We consider a three dimensional, two component, undoped, electron-hole system interacting with a trans- 
verse classical electromagnetic field along with a source for the electromagnetic field in addition to the usual 
coulombic carrier-carrier interaction. It is important to include other terms such as electron-phonon interac- 
tion as this is the dominant mechanism of energy relaxation. Impurity scattering is relatively less important 
for the ultra-pure sample that we consider. Further, it has been argued in the literature that only carrier- 
carrier scattering contributes significantly to the density dependence of the scattering rates. Intuitively, this 
is understandable as the electrons can independently emit phonons at a fixed rate irrespective of how many 
electrons are participating in the process. By contrast, the carrier-carrier scattering rate depends on how 
many carriers there are to scatter off from. We shall verify these expectations rigorously later on. 

We assume that the system is in thermodynamic equilibrium at temperature T = before time t = (we 
use natural units throughout, c = h = 1). The fields are assumed to have been turned on at t = and last 
for a duration of tx = 250/s.. The excitation energy coincides with the band gap tox = E g — Eb, where 
Eg and E^ are the band gap and exciton binding energy, respectively. Fig.l is a schematic illustration of 
the conditions under which the simulations are performed. The semiconductor Bloch equations (SBE for 
short) (see for example Ref. |33| ) are nothing but the equations of motion for the equal time component of 
the Green functions of the system. The equations in the form given below were first derived in Refs. ^ |3l| 
(Our notation is slightly different though). The Green functions in the equation below are proportional to 
the inter-band polarization and the electron and hole density distributions. 



II. SEMICONDUCTOR BLOCH EQUATIONS 



g hh (kt) = i<^(k^(ki)) = i f(k,t) 



g he (kt) = i(i> h (kt)ij e (te)) 



ij-ghhikt) = 2Re(Sl(kt)g* he (kt)) + R hh (kt) 



i—g he (kt) 



Q{kt){i - 2g hh {kt)) + (e h (k) + e c (k) - 2^{kt))g he {kt) + R he {kt) 




k' 




(1) 



k' 
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Here we have tacitly assumed that (^(ki) is equal to g ee (ld). This is an approximation that we justify 
rigorously later. The collision terms Rhh(kt) and i?^ e (kt) are in general, extremely complicated functions of 
the full Green-functions of the system. (Collisionless or Hartree-Fock(HF) approximation therefore means 
setting Rhh^t) and Rh e (kt) to zero.). Perhaps the simplest approximation for these terms would be to use 
constant relaxation rates: Rhh(kt) = —i ^^(k^/Ti, i?/j e (kt) = — i gh e (kt)/T2. The polarization dephasing 
rate is defined as 7 = I/T2. and 1/T\ is the population relaxation rate. These aprroximations are too crude 
for most purposes unless one is interested exclusively in studying coherent phenomena. The approximation 
for Rhh (kt) in terms of a constant population relaxation rate suggests that the total kinetic energy of the holes 
decays with the same rate as the total number of holes. We know that this is not true, because when phonons 
are present the kinetic energy relaxation takes place on a time scale of picoseconds whereas population 
relaxation through recombination occurs on a time scale of nanoseconds. Therefore the approximation for 
Rhh{kt) in terms of a constant relaxation rate is particularly bad. A more sophisticated approach would 
involve the use of Fermi's Golden rule to evaluate the collision terms. Here, the relaxation terms have a 
complicated dependence on the polarization and densities. However, the rates Rhhi^-t) and i?he(ki) still 
depend on the polarization and densities at time t and not on all previous times. Higher in the hierarchy of 
theoretical sophistication is the nonequilibrium Green function approach. Here, the collision terms depend on 
the full Green function of the system and therefore effectively, on the polarization and densities at all times 
previous to time t at which the collision rates are evaluated. Such a depndence of the collision terms on the 
history of the system is known as the memory effect. Each approach has its pros and cons. While the method 
of using constant relaxation rates has in its favor, simplicity, it contains phenoipenological parameters that 
have to be fitted from experiment. In fact the phenomenon of the photon echot£l has been succesfully used 
to measure the polarization dephasing rate {jr)- The use of Fermi's Golden rule works well for studying 
relaxation via phonons but is ambiguous when applied to carrier-carrier scattering (for reasons mentioned 
in the introduction and for reasons that will become clearer later.). In light of these arguments we take 
an alternative route. The collision terms in the Boltzmann equation occur as the product of the collision 
self-energy and the full Green function. The collision self-energy in turn depends on the full Green function. 
When the external fields are absent the collision terms are identically zero. Therefore, it is sensible to 
expect a power series expansion for the collision terms in powers of the external fields. In order to simplify 
matters further, we assume that for the purpose of evaluating these terms, it is sufficient to retain only 
the coupling of the carriers to the electromagnetic field, as the latter is responsible for generating all the 
dynamics. Therefore, the zeroth term in this expansion is just the collision self-energy evaluated using the 
free Green function multiplied by the free Green function which turns out to be identically zero. The next 
term would involve the first power of the total vector potential multiplied by a linear response term which 
we evaluate. A more precise quantitative meaning of this approximation will be gijven in the next section. 
The screened fock self-energy ( the GW approximation as it is sometimes called )E3 involves the evaluation 
of the screened coulomb interaction. To this end, we introduce the concept of nonequilibrium dynamical 
screening. 

A. Nonequilibrium Dynamical Screening 

The RPA is obtained by linearizing the Hartree equations of motion for the Green function in the external 
field. Alternatively, it may be obtained by replacing the full Green function by the Hartree Green function 
in the shielded potential approximation as discussed in Ref. This is exactly what we do, except that in 
our case the Hartree Green function lacks time translation invariance and therefore it is not possible to solve 
the resulting integral equation for the sheilded potential by transforming to the frequency domain. This will 
force us to seek an expansion in terms of a small parameter. 

In order to derive the relevant formulas for the screened coulomb interaction, we proceed as follows. The 
equation of motion for the mean-field (Hartree) Green function of the system, in the presence of an external 
classical electromagnetic field may be written as: 

(*A _ ff (l))G(l, 1') = 5(1 1) + £A(l).p? M G(l, 1') + ^|A(1)| 2 G(1, l'). (2) 

Here and henceforth, we use the abbreviation 1 = (xiti). This Green function will be used as the starting 
point for derivation of the formulas for the screened coulomb interaction. Note here that the Hartree-Coulomb 
terms are absent because of charge neutrality. The Green function in the above equation is defined as 
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G(l,l ) 



G ee (l,l) G* eh (l,l) 
G fce (l,l') G hh (l,l') 



and 



M 



1 

1 



(3) 



(4) 



Here the contour-ordered Green functions are defined as, 

-i{T c S ^(1)4(1')) 



G hh (l,l) 



T r S 



G ee (l,l) 



-i{T c S y>t(l)^ e (l')) 



G* eh (l,l) 



-t(T c S ^1(1)^(1')) 
T c 5 



G he (l,l ) = 



-i(T c s Mi)M±')) 

T C S 



(5) 



Here £ is the S-matrix that involves the part of the hamiltonian that generates the external field. Using the 
same methods outlined in Ref. 27 the screened coulomb interaction ug(l,3) satisfies: 



«s(l,3) = »(l-3)-i(25+l) / d2 d4 v s (l, 2)F(2, 4)v(4 - 3) 



F(2, 4) = G he (A, 2)G* eh (2, 4+) + G hh (A, 2)G hh (2, 4+) + G ee (4+, 2)G ee (2, 4) + G^(4+, 2)G, ie (2, 4). (6) 

We can now expand us (1,3) in powers of the vector potential by expanding the Green-function in powers 
of the same using the equation of motion for the mean-field Green function. The dimensionless quantity 
that appears in this expansion is denoted by A(kt). The expansion in powers of the vector potential is an 
expansion in powers of the dimensionless quantity defined below. 



A(k,t) = dt —A ext (t ) • p cv exp{-ie{k)t ), 
Jo m c 



A = — A ext .pc V Tx- (7) 

mc 

Here t is some arbitrary time after t — and A ext is the amplitude of A ext (t) and Tx = 250/s. is the pump 
duration. We require the magnitude of A(k, t) quantity to be small compared with unity for all values of 
e(k) = + Eg which is the interband transition energy at wavevector k. The expansion in powers of A is 
synonymous with working at low-densities. The leading terms in the expansion for the Green functions can 
be derived by expanding the Green function in powers of the vector potential. 

G(l, l') = G (l, l') + J d3 (^) 4(3) + \ J d3 J dA (ig^) A(3)A(4). (8) 

The derivatives are computed in the usual manner by first inverting Eq. (2) to obtain an equation for the 
inverse G~ 1 (l, 1 ). 

G-\l, 1) = (< A - Ho (l))6(l - 1) - £a(1).p? M 6(1 1) ^|A(1)| 2 ^(1 - l') (9) 
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and then by using relations such as, 

SG(l,l') = - J d3 J d4 G(l,3)(<fG _1 (3,4))G(4,l'). (10) 

Such an approach leads us to expressions for the leading terms in the expansion for the Green functions of 
the system. 

= i^0 ck (xi)(/>* k (x 1 Oexp(ie c (k)(ti -t^)), 

k 

G>(l,l') = -z^ ( /» ck (x 1 )^ k (x 1 0e^e c (k)(t 1 -t 1 0)A(kt 1 )A*(kt 1 0, 

k 

Glhi 1 ' 1 ') = -^^ ck ( Xl )^k( x i') e 2 ; P( ie c(k)ti)ea;p(ie, l (k)t 1 /)A(kti), 

k 

G lh (!) ! ) = -^^k(xi)Ck( x i') ea; P( ie c(k)ii)ea;p(ie /l (k)f 1 /)A(kt 1 '), 

k 

GfceC 1 ' 1 ') = X)^" k ( Xl ^ck( X l') ea;p ( _ie ' l ( k )* 1 ) ea; ^ _iec ( k )*l') A *( k *l')' 
k 

G he( 1 ' 1 ) = ^k(xi)0* k (x 1 - )exp(-ie h (k)t 1 )exp(-ie c (k)t l/ )A*(k*i), 

k 

= i^M^iW^Mexpi-iehimt - V))A(kt 1 0A*(kt 1 ), 

k 

Gfch(M') = -i^^k(xi)^ k (x 1 0ea:p(-ie h (k)(i 1 - i^)). (H) 

k 

From the above formulas it is obvious that the leading contribution to the collision part of the screened 
coulomb interaction depends quadratically on the vector potential. 

VS (1,2) = J d 3 ^(^||M(3)i(4) (12) 

Assuming space translational invariance (in the one conduction band and one valence band model), the 
(fourier transform of the collision part of) the screened coulomb interaction is 

< (hM) = + l)v\ ^{-Fi(k, q; h, t 2 )A(k + q; ii)A*(k; i 2 ) 



-F*(k, q; t 2 , ir)A(k; ti)A*(k + q; i 2 ) 
+F h (k, q;t 2 -ti)A(k;ti)A*(k;t 2 ) 
+ F e (k, q; i 2 - ii)A(k; ii)A*(k; i 2 ), } (13) 
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v<(t 1 ,h) = -(v>^(t 1 ,t 2 ))*, (14) 
F 1 (k,q;t 1 ,t 2 ) = exp(i(e h (k + q) - e h (k))t 2 )exp(i(e c (k + q) - e c (k))ti), 
F h (k,q;t 2 - ii) = exp(i(e h (k + q) - e ft (k))(i 2 - ii)), 

F e (k, q; i 2 - t x ) = exp(i(e c (k - q) - e c (k))(t 2 - ii)), 

w q = in two dimensions, and u q = in three dimensions. The spin S 1 = 1/2 for fermions. e c (k) = 
+ Eg and e^(k) = and the effective masses are m* = 0.067 m e and m* h = 0.5 m e and e = 12.4. 

From the above formulas for the screened coulomb interaction it is clear that for small q the screened 
coulomb interaction v>(ti,t 2 ) is proportional to in three space dimensions and independent of q in two 

space dimensions. The sum over all q's has a q 2 in the measure in 3D which cancels the divergence coming 
from 4y. In 2D the measure has a q which mutiplies a constant and vanishes for small q. Therefore, in 
both these cases we obtain a finite collision cross-section: an important and essential feature of any theory 
of screening. The above formulas differ from the naive generalisation of the RPA in some respects. Firstly, 
the traditional RPA is invoked in a slightly different context. Here, the electrons are treated as quantum 
objects immersed in a uniform positive background which does not participate in the dynamics. The uniform 
positive background ensures charge neutrality. The RPA then is just the linear response of such a system to 
weak external longitudinal (scalar) potentials. In our case, however, the external fields are weak transverse 
electromagnetic fields, but more importantly, the carriers responsible for the screening are generated out of 
the vacuum by these weak external fields themselves. So it is important to use a nonequlilibrium formalism 
from the start. The connection between the formulas derived above and RPA in the steady state limit 
becomes clear when one realises that in the steady state limit the Green functions possess time translation 
invariance and therefore Eq. (6) can be fourier transformed and inverted algebraically to yield an expression 
for the screened coulomb interaction. This screened interaction is expressed in terms of the full (steady 
state) Green function of the system. (Physically, the steady state limit can be realised by applying an 
infinitely weak field for an infinitely long time.) When the full Green function is replaced with the the 
mean-field Hartree Green function, we recover the familiar RPA. If we assume further that the densities 
involved are small, it is easy to see that the collision part of the screened coulomb interaction is quadratic 
in the unscreened potential. This is equivalent to iterating the integral equation (Eq. (6)) once. The steady 
state limits of Eqs. (13) and (14) are therefore identical to the low density limit of RPA. However, the RPA 
(i.e., the replacement of the full Green function by the Hartree Green function in Eq. (6)) is known to be 
exact only in the ultra high density limit and its low density limit is not likely to capture all the physics. In 
light of the complexity of the problem we are dealing with here, and all the other approximations that we 
have made, it is reasonable to proceed further, after pointing out these pitfalls to the reader. 

Quasistatic screening is the popular choice in the literatureQ. But it is well known that this type of 
screening grossly overestimates the importance of screening yielding an exponentially decaying Yukawa type 
of potential in real space. The naive generalisation of the RPA involves replacing in the formula for the 
equilibrium RPA, equilibrium densities(and polarizations) by nonequilibrium ones. When this is done, it 
not not clear what the meaning of the frequency that appears in the dielectric function is. Nor is it clear 
how that screened coulomb interaction should be coupled with Fermi's Golden Rule. Our approach is more 



systematic and perhaps more appealing. The approach by authors in Ref. 36 comes closest to the one just 
discussed above. The authors however do not solve the SBE including the detailed nonequlibrium screening 
that they introduce. 

B. Collision Terms due to Screening of the Coulomb Interaction 

The collision terms are evaluated using the same scheme that was used in the evaluation of the screened 
coulomb interaction. Rh e (kt) is expanded in powers of A(kt). The collision self-energy is E c (l,l ) = 
w s (l,l')G(l,l') Defining 



(i 



*»'>-(&<M : ) S'a.l')) =/«^*Wl')-/^G(l,2)E e (2,l'), 
the relevant component of i2(l, 1) namely, 

R he (x 1 t- L ,x l/ t 1 +) = ^^„ k (xi)^* k (x 1 /)i?/ le (kti) 

k 

yields the necessary collision term. An expansion of R(l, 1 ) in powers of A implies that the collision term 
Rhe(kt) is proportional to A 3 . 




U 1 (k,q;t,t 2 ) = v>(t,t 2 )exp(-i(e h (k - q) + e c (k))t)exp(-i(e c (k - q) - e c (k))i 2 ) 

+«q (i2,i)exp(-i(e h (k - q) - e h (k))t 2 )exp(-i(e c (k - q) + e/»(k))i), 
C/ 2 (k,q;t, t 2 ) = «q (i, t 2 )exp(i(e/ l (k - q) - e h (k))t 2 )exp(-i(e h (k - q) + e c (k))i) 

+ «<(t 2) t)esp(-i(c c (k) - e c (k - q))i 2 )ea*(-i(e;,(k) + e c (k - q))t), (15) 

and Rhh(k, t) and i? ee (k, t) are of the order of A 4 and are therefore neglected as far as coulomb scattering is 
concerend. 

C. Collision Terms due to Phonons 

For treating the electron-phonon collisions we must expand both Rhh{kt), R ee (kt) and Rhe(kt) in powers 
of A(k£). The expansion yields the following from for the collision terms. 

R he {kt) = - f dt 2 1 V|M q | 2 sin(^ LO (i 2 -t))[ea;p(z(e c (k-q)-e c (k))i 2 ) 
Jo v q 

exp{-i{e c (k - q) + e/,(k))t)A*(kt 2 ) 
-earp^e/^k) - e h (k - q))i 2 ) 
«cp(-i(e c (k - q) + e h (k))i)A*(k - qt 2 )] 

_ / dt2 77 y]l^q| 2 sin(u; LO (i-i 2 ))[ea;p(i(e ?1 (k-q) - e h (k))t 2 ) 
Jo V q 

ea;p(-i(e ft (k - q) + e c (k))i)A(k - qi 2 ) 
-ea:p(i(e c (k) - e c (k - q))i 2 ) 
e:cp(-i(e c (k) + e h (k - q))i)A(ki 2 )]A*(k - qt)A'(ki) 
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and similarly for RhhQct). 

R hh (kt) = -y^M^ 2 sin(u LO (t - t 2 ))[-exp(i{e c {k) -e c (k-q)i 2 ) 
q 

erp(i(eh(k) - e h (k - q))i)A(kf 2 )A*(k - qi) 
+exp(i(eh(k - q) - e h (k))i 2 ) 
ea;p(i(efc(k) - e^k - q))t)A(k - q< 2 )A*(k - qi)] 
+c.c. 

It can be seen that the leading contribution Rh e due to phonons is proportional to A and the next term is 
proportional to A 3 . In the term proportional to A it appears that we have to use the total vector potential 
and not the external vector potential. This is because corrections to A resulting from such a replcement 
is proportional to A 2 which is larger in magnitude than the next leding term in Rh e which is proportional 
to A 3 . But in fact, this apparent inconsistency is not so severe, as we have found that for times when the 
exciting fields are present, the two A's are within 1% of each other. However, for times when the external 
fields are absent, it is important as a matter of principle to include this correction. However, we neglect 
this correction because of a more basic reason, namely when the external fields are off the fluctuations in 
the vector potential are going to be comparable to the expectation value of the same. Since we are anyway 
going to treat the fields classically, we are obligated to ignore the difference between the two A 's as well, 
after the external fields are turned off. 

The fact the the electron and hole effective masses are different means that Rhh(kt) and i? ee (ki) are not 
exactly equal. This is inconsistent with charge conservation. A remedy is to use the reduced mass instead of 
the effective masses thereby forcing charge conservation. This approximation will not affect the qualitative 
results and has the added appeal of obeying a conservation law. 

D. Scaling Relation of the Dephasing Rate with Density at Ultra-Low Densities 

The theory of the semiconductor photon echo is fairly well develope The central feature of the theory 
is that the energy of the echo signal is known to vary exponentially with the the time delay r between the 
two pulses invloved in the echdla. Quantitatively, E(t) a exp(—T/T ec h ) and the echo time is related to the 
polarization dephasing rate 7 through T ec ho — j*-. The echo time can be measured by measuring the echo 
energy as a function of the temporal separation of the pulses. The polarization dephasing rate 7 is a phe- 
nomcnological constant that accounts for the decay of the polarization due to various relaxation processes. 
From our analysis it is possible to predict with accuracy the scaling relation of this phenomcnological dephas- 
ing rate with the total carrier density at ultra-low densities. We have argued that the density dependence of 
the relaxation rate is solely due to carrier-carrier scattering. Our analysis shows that i?; le (ki) is proportional 
to A 3 for coulomb scattering and proportional to A for phonon scattering. Further, the total density which 
is related to <7wi(ki) scales as A 2 and gh e (kt) is proportional to A for small A. Therefore the polarization 
dephasing rate at ultra-low densities is proportional to the density (from carrier-carrier scattering) and also 
has a added part independent of the density due to phonons. A critical evaluation (and a measure of the 
smallness of the density) of this conclusion as well as a comparison with the results of other authors will be 
given at the end. 

III. COMPUTATIONAL PROCEDURE 

Having obtained the collision terms, we are now ready to solve the semiconductor Bloch equations nu- 
merically. We follow closely the approach outlined in Ref. |3^. We assume an isotropic distribution and 
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perform all angular integrations analytically. The singular coulomb terms in the HF part are particularly 
troublesome and need special care in their evaluation. This is a well-known approach outlined in Ref. |33j. 
For example, a sum of the type shown below can be rewritten as described: 

E «k-k'/(*')= E «k-w(/(*') -/(*))+/(*) E u k- k <- 

k^k' k/k' k/k' 

The last part can be evaluated anaytically and in the first part the angular intagrations performed analyti- 
cally. We also assume a cut-off for the wavevector k max . The above sum may then be written as 



E 

k^k' 



'k-k 



/(*) 



(Aire 2 ) 
(2tt)* 



dk (—)ln( 
k 



(k + k) 



(k' - k) 



)(/(*)-/(*)) 



(47re 



r/(*)(^)^(*S»ax-* 2 )M 



) -|- k k n 



The cutoff k max is chosen arbitrarily and its validity is justified aposteriori. For example, for the case of 



excitation at the excitonic absorption edge we choose k„ 



This choice is justified on the grounds 



that the tail of the occupation probability distribution becomes close to zero as the wavevctor approaches 



We choose A ext (t).p V c = A pump p vc cos{oj pump t)9{Tx - t) and A p 



7Z £— v. Here, p vc is the 



• J pump L J L ' V ' A v J oia.lv* ^ipump — (p uf ^ e Tx ) ' 

inter-band momentum matrix element, tx — 250/s. is the pump duration. We can vary the intensity by 
tuning the dimcnsionlcss parameter Ao. For the numerical solution of the semiconductor Bloch equations we 
use the standard Runge-Kutta procedure. The screened coulomb interaction is first numerically evaluated 
and tabulated in an array before proceeding to the Runge-Kutta algorithm. It is important to use the 
rotating wave approximation described in detail in the seminal work of R. Binder et. al.Lj in order to 
get sensible results. The collision terms are partly evaluated analytically, tabulated and then used in the 
computations. 



IV. RESULTS AND DISCUSSIONS 



The main results of our formalism are shown in Figs. 2 and 3. The phenomenological dephasing rate 
7 = I/T2 is plotted as a function of the total excited carrier density. The simulations as depicted in Fig. 2a 
confirm our claim that at ultra-low density, the dephasing rate scales linearly with the density. The scale for 
the smallness of the density is set by the value of Ao- In fact, our simulation shows (Fig. 2b) that the total 
excited carrier density at ultra-low densities is given by po = 7.8Aq (in units of I0 15 cto~ 3 ). Therefore, for a 
choice of Ao = 0.1 which is small compared to unity, we find that the upper limit on the densities that we 
can work with is 7.8 x I0 13 cm -3 At these low densities, unfortunately, experimental data are hard to come 
by, because the signal becomes very weak at low densities. It is worthwhile to stress the following point. 
The fact that we are able to extract the density dependence of the dephasing rate does not mean that the 
actual relaxation via carrier-carrier scattering is simply an exponential relaxation. The quantity T2 is merely 
a useful phenomenological parameter that gives us an estimate of the rate of relaxation. It also happens to 
be directly measureable experimentaly, via four-wave mixing. Strictly speaking, we should have simulated 
the four-wave mixing experiment including the realistic Coulomb scattering, plotted the energy of the echo 
signal as a function of the time delay between the pulses, and then extracted the dephasing rate from the 
plot. This is much too difficult to carry out and therefore we have chosen the simpler route outlined here. 

It is possible to go beyond these densities by making a slight modification to our formalism. The phe- 
nomenological collision term i?^ e contains terms proportional to the powers of the polarization via the 
self-energy and the Green function. The polarization in turn decays with a time constant of T 2 . So, Rhe 
itself should decay at the same rate or faster. In our nonequilibrium formalism however, this does not hap- 
pen. The decay has to come from the polarization which is proportional to A. This A as we have defined it, 
does not decay with time. However, we know that the polarization does. So we define a new A as follows. 

A'(ki) = A(kt)exp(-t/T 2 ) 
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with a T2 to be determined later. Therefore, we replace all the A(kt) 's in the collision terms by A (kt)'s. This 
means that the actual polarization (for small Ao) is proportional not to A(k<) but to A (kt). This ensures 
that the collision terms decay with time. The time constant T 2 has to be determined self-consistently. We 
vary the value of T2 and plot the total density as a function of time(not shown). For time scales that we 
consider, ( less than 1 ps. ) no recombination occurs and so we should see a constant total density when 
the external fields are switched off. For Ao close to unity and large enough T2, we see that the theory is 
asymptotically unstable. In other words, the total density increases without saturation. This unphysical 
result comes about because of a lack of self-consistency in the formalism. In order to derive expressions 
for the screened coulomb interaction and the collision terms we have ignored the effects of collisions. This 
inconsistency manifests itself as asymptotic instability. In other words, the theory is unreliable at long time 
scales. The device of introducing T2 is precicesly to partially restore this self-consistency. The value of T2 is 
tuned from a large value to smaller and smaller values till we see that the total excited carrier density remains 
constant after the exciting field has been switched off. Any smaller value of the T2 than this critical value will 
also yield the same result, of course. However, too small a value of T2 will spoil the small time behaviour of 
our formalism which is modeled exactly. Therefore, the critical value for T2 provides the dephasing rate. The 
Fig. 3a shows the result of the simulations. The experimental data from Ref. [l6| shows that at densities larger 
than 10 16 cm -3 the polarization dephasing rate scales as dpt. 01 Pq' 33 ■ This sublinearity in the scaling is 

confirmed by our simulations. The difference between the predicted exponent and the observed is probably 
atrributable to lack of full self-consistency in the formalism and the use of a simple two-band model. The 
sublinearity originates because of a saturation effect on the dephasing rate tjt brought about by imposing 
partial self-consistency to the formalism. We have found that the dephasing rate at relatively high densities 
(> 1 x 10 15 cto -3 ) is given by 

1 -imps- 1 ) = 4.0(p ) - 48 



1 2 

where po is in 10 15 cm~ 3 . In addition, there is a additive contribution of 5 ps -1 to -jr. independent of 
the density (for densities in the range 10 12 to 5 x 10 15 cm -3 ) coming from phonons (not shown in the 
figs.). Therefore, the carrier-carrier dephasing rate becomes equal to the carrier-phonon dephasing rate at a 
density of 2 x 10 15 cm~ 3 . These results are also confirmed by experiments in Refs. H and The conventional 
wisdom that carrier-carrier scattering is much faster than carrier-phonon scattering is true at densities much 
higher than those considered here. The dependence of the density on the parameter Ao is shown in Fig. 3b. 
The dependence is non-parabolic as Ao approaches and exceeds unity. The dephasing rate for phonons 
is computed in an manner analogous to the method used for carrier-carrier scattering. The fact that the 
dephasing rate is independent of the density in the density range considered here suggests that the term 
that is proportional to A 3 contributes insignificantly. We have confirmed this by plotting various quantities 
with and without this term and registering only a marginal change. The explanation is, a possible accidental 
cancellation of the various components of this term. In fig. 4 we see that for a pulse Ao = it, in the absence 
of collisions (shown in dotted lines) the phenomenon of quantum beats which is the oscillation of the hole 
occupation probability in time after the external field has been turned off ( tx = 250/s.). This is due to the 
superposition of discrete excitonic frequencies resulting in the beating behaviour. The presence of carrier- 
carrier scattering (in solid) results in the decay of the quantum beats with a rate given by ^- = 21.9/ps.. 
The value of Ao = tt is much too large compared to unity and significantly beyond the range of applicability 
of our formalism, and therefore the value of T 2 quoted above should not be taken too literally. The large 
value of Ao ensures a large value of y- and therefore emphatically illustrates the phenomenon of dephasing. 



V. CONCLUSIONS 



We have been able to construct an elaborate theory of screening in the presense of external classical time 
varying electromagnetic fields. We have obtained explicit formulas for the screened coulomb interaction, 
and demonstrated that they produce finite collision cross-sections in both two as well as three dimensions. 
We have also explained the range of validity of these formulas. Using the theory we have demonstrated 
the following: The density dependence of the polarization dephasing rate is due to carrier-carrier scattering 
alone, in the range of densities considered, and the energy relaxation is due mainly to phonons. Phonons 
contribute to the dephasing rate as an additive term independent of the density. We have also derived scaling 
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relations of the polarization dephasing rates with density. We have also provided quantitative estimates of 
the dephasing rates and demonstrated that they agree qualitatively with experiments. In particular, we have 
demonsrtated that the dephasing times at low densities due to carrier-carrier scattering is in picoseconds 
and not femtoseconds as is sometimes believed. We have also studied the time evolution of the probability 
distribution and demonstrated the phenomenon of dephasing of excitonic quantum beats as a result of 
carrier-carrier scattering. Further imporovements should involve the use of more realistic band structures 
such as the inclusion of valence subbands. 
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Fig.l Figure demontrating the conditions under which the simulations are perfeormed. 

Fig. 2a Linear scaling of the dephasing rate with density at ultra-low density. 

Fig. 2b Parabolic dependence of the the excited carrier density on Aq. 

Fig. 3a Sublinear scaling of the dephasing rate with density at higher densities. 

Fig. 3b Non-parabolic dependence of the the excited carrier density on Ao for higher densities. 

Fig. 4 Dephasing of the excitonic quantum beats as a result of carrier-carrier scattering. 
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